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Abstract 

Relying on recent advances in the theory of entropy solutions for nonlinear (strongly) 
degenerate parabolic equations, we present a direct proof of an L 1 error estimate for 
viscous approximate solutions of the initial value problem for 

d t w + div(V (x)f(w)) = AA(w), 

where V = V(x) is a vector field, / = f(u) is a scalar function, and A'(-) > 0. The 
viscous approximate solutions are weak solutions of the initial value problem for the 
uniformly parabolic equation 

d t w £ + div(V(x)f{w e )) = A(A{w e )+ew e ) 1 e > 0. 

The error estimate is of order y/s. 

1 Introduction 

In this paper we are interested in certain "viscous" approximations of entropy solutions 
of the initial value problem 

d t w + div(V(a?)/(u7)) = AA(w), (x, t) e Qt, 

w(x,0) = w (x), x E M. d , (1.1) 

where Q T = R d x (0, T) with T > fixed, u : Q T -> R is the sough function, V : R d -> R is 
a (not necessarily divergence free) velocity field, / : M — > R is the convective flux function, 
and A : R — > R is the "diffusion" function. For the diffusion function the basic assumption 
is that A(-) is nonincreasing. This condition implies that QLlj) is a (strongly) degenerate 
parabolic problem. For example, the hyperbolic equation dtw + divf V(x)f(w)) = is 
a special case of ()1.1|) . Problems such as (|1.1[1 occur in several important applications. We 
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mention here only two examples: flow in porous media (see, e.g., |Sj) and sedimentation- 
consolidation processes 

Since A(-) is merely nondecreasing, solutions are not necessarily smooth and weak solu- 
tions must be sought. Moreover, as is well-known in the theory of hyperbolic conservation 
laws, weak solutions are not uniquely determined by their initial data. To have a well- 
posed problem we need to consider entropy solutions, i.e., weak solutions that satisfy 
a Kruzkov- Vol'pert type entropy condition. A precise statement is given in Section |2] (see 
Definition [J). For purely hyperbolic equations this entropy condition was introduced by 
Kruzkov ^S] and Vol'pert For degenerate parabolic equations, it was introduced by 
Vol'pert and Hudjaev [22j. 

Following Carrillo 5 , Karlsen and Risebro proved that the entropy solution of 
(as well as a more general equation) is unique. Moreover in the L°° (0, T; BV (R d )) 
class of entropy solutions, they proved an L 1 contraction principle. Existence of an 
L°°(0,T;BV(M. d )) entropy solution of follows from the results in Vol'pert and Hu- 
djaev |22| or Karlsen and Risebro |12j (the latter deals with convergence of finite difference 
methods) . The proof in of uniqueness and stability is based on the "doubling of vari- 
ables" strategy introduced in Carrillo [H] (see also Chen and DiBenedetto [0]), which in 
turn is a generalization of the pioneering work by Kruzkov |15j on hyperbolic equations. 
Related papers dealing with the "doubling of variables" device for degenerate parabolic 
equations include, among others, Carrillo 4 , Otto ^F3\, Rouvre and Gagneux [201, Cock- 
burn and Gripenberg [Jj Burger, Evje and Karlsen PJE], Ohlberger JH]; Mascia, Por- 
retta, and Terracina ^7j> Eymard, Gallouet, Herbin and Michel and Karlsen and 
Ohlberger [T3] . 

In this paper we are interested in certain approximate solutions of (jl.lj) coming from 
solving the uniformly parabolic problem 

d t w £ + dw(V(x)f(w £ )) = AA £ (w £ ), (x, t) £ Q T , 

w £ (x,0) = w (x), x£R d , (1.2) 

where A £ (w £ ) = A(w £ ) + ew £ , e > 0. We refer to w £ as a viscous approximate solution 
of 1)1 -lj) . Convergence of w £ to the unique entropy solution w of (|1,1|) as e 1 follows from 
the results in Vol'pert and Hudjaev [22]. Our main interest here is to give an explicit rate 
of convergence for w £ as e [ 0, i.e., an L 1 error estimate for viscous approximate solutions. 

There are several ways to prove such an error estimate. One way is to view it as 
a consequence of a continuous dependence estimate. Combining the ideas in ^3] with 
those in Cockburn and Gripenberg [Jj, who used a variant of Kruzkov's "doubling of 
variables" device for ()1.1|) with V = 1, Evje, Karlsen and Risebro 9 established an explicit 
"continuous dependence on the nonlinearities" estimate for entropy solutions of (jl.lj) . 
A direct consequence of this estimate is the error bound 1 1 w £ — w \ \ L i (q t ) = O , at least 
when w £ ,w belong to L°° (0, T\ BV fR")) and V is sufficiently regular. Unfortunately the 
techniques employed in |S] require that one works with (smooth) viscous approximations 
of (jl.l|) . The proof in [0] (as well as the one in jjj) did not exploit the entropy solution 
"machinery" developed by Carrillo [3]. 

The main purpose of this work is to show that one can indeed use the "doubling of 
variables" device to compare directly the entropy solution w of (|1.1|) against the viscous 
approximation w £ of (|1.2|) . Hence there is no need to work with approximate solutions 
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of Although our proof is of independent interest, it may also shed some light on 

how to obtain error estimates for numerical methods. Most numerical methods (related 
to this class of equations) have (|1,2[) as a "model" problem and, in this context, the 
size of e designates the amount of "diffusion" present in the numerical method. A step 
in the direction of obtaining error estimates for numerical methods has been taken by 
Ohlberger |18j with his a posteriori error estimate for a finite volume method. We will 
in future work use the ideas devised herein to derive error estimates a priori for finite 
difference methods. 

The rest of this paper is organized as follows: In Section |2] we state the definition of an 
entropy solution and the main result (Theorem^). Section|31is devoted to the derivation of 
certain entropy inequalities for the exact entropy solution and its viscous approximation. 
Equipped with these entropy inequalities, we prove the error estimate (Theorem ^) in 
Section 0J 

2 Statement of result 

Following EH we start by stating sufficient conditions on V = (Vi, . . . , V^), /, A and 
uo to ensure the existence of a unique L°°(0, T; BV(R d )) entropy solution of 1)1. ljl : 

V e (L°°(R d )) d n (Lip(M d )) d ; divV G BV(R d ); 
f G Lip loc (M); /(0) = 0; 

A G Lipi oc (R) and A(-) is nondecreasing with ^4(0) = 0; 

u G L°°(R d ) n BV(R d ). (2.1) 
Note that the first condition in ()2.1|) implies 

In (|2.1j) and elsewhere in this paper the space BV(R d ) is defined as 
BV{R d ) = [g G L l {R d ) : \g\ Bv{Rd) < oo} , 

where M^y^^ denotes the total variation of <?, i.e., g G BV(R d ) if and only if g G L 1 (R rf ) 

and the first order distributional derivatives of g are represented by finite measures on R d . 
Equipped with (|2.1j) we can state the following definition of an entropy solution: 

Definition 1 (Entropy Solution). A function w(x,t) is called an entropy solution 
of JIIJ if 

(i) w G ^(Qt) PI L°°(Qt) PI C (0, T; L 1 (R d )), 

(ii) A(w) G L 2 {Q,T-H l (R d )), 

(iii) w(x, t) satisfies the entropy inequality 

w - k\d t 4> + sgn(w - k) [V(x)(f(w) - f{k)) - VA(w)] • V0 

Qt 

- sgn(w - k)divV(x)f(k)4>) dtdx>0, Vfcel, (2.2) 
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for all nonnegative (ft £ C^(Qt) and 
(iv) ||u>(-,t) — WolliifiKrf) as i J. (essentially). 

Note that, if we take k > esssup w(x,t) and k < essinf w(x, t) in (|2.2j) . then an ap- 
proximation argument reveals that 

wfa + [V(x)f(w) - VA(w)] ■ V^) dtdx = (2.3) 

holds for all (j) £ H 1 {Qt)- Let (•,•) denote the usual pairing between i/ _1 (lR a! ) and 
if 1 (M d ) . From ((231) we conclude that 

d t w e L 2 (Q,T;H- 1 (R d )), 

so that 

-J (d t w,<p)dt+Jj([V(x)f(w)-VA(w)]-V(j?Jdtdx = 0, V € H 1 {Q T ). (2.4) 

Qt 

In other words an entropy solution iu(x, i) of is also a u>ea& solution of the same 

problem. 

In this paper we are interested in comparing the entropy solution w of against 
the weak solution w £ of the viscous problem 1)1. 2JI . From the results in Karlsen and Rise- 
bro 12 or Vol'pert and Hudjaev [221 there exists a weak solution w £ £ L°°(0, T; BV(R d )) 
of Q1.2[) . Since j4 £ (-) is increasing, the uniqueness result in Karlsen and Risebro ^3] (see 
also Remark n herein) tells us that this weak solution is in fact a unique solution. More- 
over from the energy estimate we conclude that w £ G £ 2 (0, T; H 1 (IR^)) . Of course, if V, 
f, A, uq are smooth enough, one can prove that the weak solution w £ of (|1.2jl is actually 
a classical (C 2,1 ) solution. See, e.g., Vol'pert and Hudjaev [22] • Here it will be sufficient 
to know that w £ belongs to L 2 (0,T; H 1 (R d ) ) (not C 2 ' 1 ). 

We are now ready to state our main theorem: 

Theorem 1 (Error Estimate). Suppose that the conditions in (|2.1[) hold. Let w £ 
L°°(0, T; BV (R d ) ) be the unique entropy solution of ([□} and let w £ £ L 2 (0,T; H 1 (R d ) ) n 
L°° (0, T; 2?V (R rf )) 6e i/ie unique weak solution of (|1.2|) . T/ien i/tere exists a constant C, 
independent of e, such that 

W -A\v(Q T )<CV~e- (2.5) 
3 Entropy inequalities 

In Section 0] we follow the uniqueness proof of Carrillo 5 j to obtain an estimate of the 
difference between w £ and w. To this end it will be necessary to derive two entropy 
inequalities for the exact solution w and two approximate entropy inequalities for the 
viscous solution w £ . The purpose of this section is to derive these inequalities. (See 
Lemma |21 and Lemma 01 below.) 
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Note that, differently from the pure hyperbolic case JH]) we need to operate with 
one additional entropy inequality (actually an equality for the exact solution w) taking 
into account the parabolic (dissipation) mechanism in the equation. Hence we introduce 
a set H corresponding to the regions where A(-) is "flat" and (jl.lj) behaves hyperbolically. 
More precisely, let A^ 1 : R — > K denote the unique left-continuous function which satisfies 
A^iAiu)) = u for all u£R. Then we define 



H = |r G M : A is discontinuous at r\. 



Since A(-) is a monotonic function, H is at most countable. The dissipation mechanism in 
the equation is effective only in the (x,t) region corresponding to the complement of H. 
To prove Lemma 121 and Lemma 01 below we need the following "weak" chain rule: 

Lemma 1. Let u : Qt — ► R be a measurable function satisfying the four conditions 

(1) u g l 1 (Qt) n l°°(Qt) n c(^,T-L l {w 1 )), 

(2) n(0, •) = u G L 1 (R d ) n L°° (R d ) , 

(3) S^u€L 2 (0,r;fr- 1 (R d )) and 

(4) A{u) £ L 2 (0,T;H 1 (R d )). 

For every nonnegative and compactly supported <p £ C°°(Qt) with (p\t=o = 4>\t=T = we 
have 



(d t u^{A{u))(^dt = jJ(J ^{A{i))dxi\4>tdtdx, ke 



Qi 

where ip :R -^R is a nondecreasing and Lipschitz continuous function. 

The proof of Lemma^is very similar to the proof of the "weak chain" rule in Carrillo jS] 
and it is therefore omitted. See instead |13j . 

The following lemma, which deals with entropy inequalities for the exact entropy solu- 
tion w, is a direct consequence of the very definition of an entropy solution. 

Lemma 2. The unique entropy solution w of ([l.ljl satisfies: 

(i) For all k G IR and all nonnegative (j) G C^°(Qt) we have 

E h ^(w,k,(/)) > 0, (3.1) 

where 

E h y p (w, k, <£)■= jj (\w - k\d t <j> + sgn{w - k) [V(x)(f(w] - f{k)) 
Qt 

- VA(w)} ■ - sgn(u> - k)divV(x)f(k)(j)\ dt dx. (3.2) 
We refer to (JSHJ) as a hyperbolic entropy inequality. 
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(ii) For all k such that A(k) ^ H and all nonnegative <p S C^°(Qt) we have 

E par (w,k,0) =0, (3.3) 

where 

E^(w, k, <j>):= ff (\w - k\dt<j> + sgn(u; - k) [V(x)(f(w) - f(k)) 



- VA(w)] ■ V0 - sgn(w - k)divV(x)f(k)<f>) dtdx 

-lim J J \V A(w)\ 2 sgn'^ (A(w) - A(k))(j) dtdx. (3.4) 
Qt 

In (|3.4|) ( and elsewhere in this paper) sgn^ is the approximate sign function defined 
by 

Jsgn(r) if\r\ > rj, 
sgn (r) ■=< . .... rj > 0. (3.5) 

[T/77 if \t\ < rj, 

We refer to (|3.3() as a parabolic entropy inequality. 

Proof. The first inequality (|3.1|) is nothing but the entropy condition for the entropy so- 
lution w. So there is nothing to prove. We turn to the proof of the second inequality (|3.3|) . 
which borrows a lot from Carrillo |2] (see also ^S]). In what follows we always let k and <ft 
be as in the lemma and the approximate sign function sgn^(-) is always the one defined 
in (|3~5ll . 

Since w satisfies (HJ) and [sgn^H - A(k))<f>] € L 2 (0,T; H 1 (R d )) , we have 
-J (d t w,sgn v (A(w) - A(k))(f^dt 

+ ff([V(x)(fH-f(k))-VA(w)]-V[sgii v (A(w)-A(k))<f>] 



Qt 



dwV(x)f(k)[sgn v (A(w) - A(k))cf\) dtdx = 0. 



Introduce the function ipr]( z ) = sgn^(z — A{k)) and note that LemmaQcan be applied so 
that 

'd t w,sgn v (A(w) -A(k))(f)^dt = j j ( j sgn r? (A(^) - A(k)) dA d t (f> dt dx . 



Hence 



J j (^£sgn v (A(0 - A{k))dAjd t (j>dtdx 
Qt 

+ JJ([V(x)(f( W )-f(k))-VA(w)] ■V[ S gn v (A( W )-A(k)) ( f>] 

Qt 

- sgn^(i4(io) - A(k))divV(x)f(k)4?j dt dx = 0. (3.6) 
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Note that since A(r) > A(k) if and only if r > k (here we make use of the assumption 
that k G "parabolic region", i.e., A(k) ^ H), sgn (A(r) — A(k)) — ► 1 as 77 J. for any 
r > k. Similarly for r < k. Consequently, as 77 [ 0, J^sgn — A(k)) d£ — > \w — k\ 

a.e. in Qt- Moreover we have | f]f sgn r? (yl(^) — A(k)) d£\ < \w — c\ G L\ oc (Qt) so that by 
Lebesgue's dominated convergence theorem 

limjj^ sgn^A^) - A{k))di S jd t ^)dtdx = J j \w - k\d t 4> dt dx . 
Qt Qt 

Next we have 

lim J j [V(x)(f(w) - f(k)) - VA(w)] ■ V[sgn v (A(w) - A{k))cp] dtdx 
Qt 

= lim jj [V(x)(f(w) - f(k)) - VA(w)] ■ Vsgn v (A(w) - A{k))<j)dtdx 
Qt 

+ lim jj [V(x)(f(w) - f(k)) - VA(w)] ■ sgn v (A(w) - A{k))V <j) dt dx 
Qt 

= lim jj V(x){f(w) - f(k))sgn' v (A(w) - A(k))VA(w)<t>dtdx 



h 



- lim jj \VA{w)\ 2 sg-a' 11 {A{w) - A{k))<j)dtdx 
Qt 

+ lim jj sgn^Aiw) - A(k))[V(x)(f(w) - f(k)) -VA(w)] -V<j)dtdx. 



Note that I\ can be rewritten as l\ = lim JJ V(x)dwQ v (A(w))(j)dt dx, where 
Q v (z) := j\gn' v (r - A(k)) (/ (A' 1 (r)) - /(A" 1 ^)))) dr 

Surely Q^fz) tends to zero as 77 j for all z G Range(A). By invoking Lebesgue's domi- 
nated convergence theorem, we conclude after an integration by parts that 



h = - lim f f ( QJA(w))V{x) ■ V4> + QJA(w))dwV(x)<f>) dt dx = 0. 
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Using that sgn(w — k) = sgn(A(w) — A(k)) a.e. in Qt (since A(k) £ H) we have 



/ 2 =lim J J sgn v {A(w) - A(k))[V{x){f{w) - f(k)) -VA(w)] -V(f>dtdx 
Qt 

sgn(w - k) [V(x)(f(w) - f(k)) - VA(w)] • V(j)dt dx. 



For the same reason we have that 



lim / / sgn n {A(w) — A{k))d\vV (x) f (k)4> dt dx = J J sgn{w — k)d\vV{x)f{k)(j)dtdx. 
Qt Qt 
Consequently, letting r] I in (|3.6|) . we obtain ()3.3|) , ■ 

Remark 1. Observe that, if A(-) is increasing, then a weak solution is automatically an 
entropy solution and hence it is unique. 

The next lemma, which deals with approximate entropy inequalities for the viscous 
solution w £ , is a direct consequence of the definition of a weak solution of (|1.2j) . 

Lemma 3. Let E hyp and F par be defined in (|3.2|) and (|3.4|) respectively. Furthermore 
define 

R visc :=e \\7w £ ■ V<p\ dt dx. (3.7) 



The unique weak solution w £ G L 2 (0, T; H 1 (R d )) n L°° (0, T; BV (R d ) ) o/ (JOJ) satisfies: 

(i) For aZZ G R and all nonnegative <f> G C^{Qt) we have 

E h ^(w e ,k,cP)>-R visc . (3.8) 
We refer to (|3.8|) as an approximate hyperbolic entropy inequality. 

(ii) For a// G K suc/t i/iai -A(A;) ^ iif and all nonnegative <j) G C^{Qt) we have 

E^(w £ ,k,<P)>-R visc . (3.9) 
We refer to (|3.9|) as an approximate parabolic entropy inequality. 

Proof. In what follows we always let k and 4> be as indicated by the lemma. The proof 
of the inequality (|3,8|) follows the proof of ()3.1|) rather closely. Since w e is a weak solution 
and [sgn v {w £ - k)(j>] belongs to L 2 (0,T; H 1 (R d )) , we have 

T 

(dtw £ , sgn v (w £ — k)<j)J dt 

+ U ([V(x)(f(w £ ) - f(k)) - VA £ (w £ )] ■ V[sgn v (w £ - k)<j>] 

Qt 

- d\wV{x)f(k)[sgn v {w £ - k)<p]^ dtdx = 0. 
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By the chain rule we obviously have 
-J (d t w £ ,sgn v (w £ -k)<pjdt 

= JJ Q sgn^(£ - k) <%J d t <t> dt dx ^ j I \w £ -k\d t <f>dtdx 

Qt 

so that 

JJ \w £ — k\dt4> dt dx 

Qt 

+ lim J j' (\V{x){f{w £ ) - f(k)) - VA £ (w £ )] • V[sgn> £ - k)<f>] 



IT 



- sgn v (w £ - k)divV{x)f(k)(f^ dt dx = 0. (3.10) 
Firstly we have 

lim JJ sgn^w/ — k) divV {x) f (k)4> dt dx = JJ sgn(u> £ — k)divV(x)f(k)4> dt dx. 



Qt Qi 
Next we have 



lim J j [V(x)( f(w £ ) - f(k)) - VA £ (w £ )] ■ V[sgn v (w £ - k)<f>] dtdx 
Qt 

= lim ff [V(x)(f(w £ ) - f(k)) - VA £ (w £ )] ■ Vsgn (w £ - k)<f>dtdx 
Qt 

+ lim JJ [V(x)(f(w £ ) - f(k)) - VA £ (w £ )] ■ sgn T; (u> £ - k)V<j)dtdx 
Qt 

= lim ff V(x)(f(w £ ) - f(k))sgn'(w £ - k)VA £ (w £ )(pdtdx 
Qt 

v ' 

h 

~ 1 ™JJ (A £ )'(w £ )\Vw £ \ 2 sgn f v (w £ - k)<j)dtdx 
Qt 

+ JJ sgn(w £ - k)[V(x)(f(w £ ) - f(k)) - VA £ (w £ )] -V^dtdx. 

Qt 

Note that I\ can be rewritten as l\ = lim ff V(x)dwQ v (w £ )(j) dt dx, where 

^Qt 

Q v (z) := [ Z S gn' v (r-k){f(r)-f(k))dr 
Jo 

= - / (f(r) - f(k)) dr - as r, j 0. 

V J min(z,fc— r/) 
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From Lebesgue's dominated convergence theorem we conclude that 

h = - lim // (q v (A{w £ ))V{x) • V0 + Q r ,{A(w £ ))divV{x)<p) dt dx = 0. 

r ll°JjQ T ^ ' 
In conclusion we have 

vf - k\4> t + sgn(u; e - k) [V(x)(f(w £ ) - f(k)) - VA £ (w £ )] ■ V(f> 



Q 



T 



— sgn(w £ — k) divV (x) f (k)4>) dtdx 



lim / / (A £ y(w £ )\Vw £ \ sgn v {w £ - k)<j)dtdx > (3.11) 



Q 



T 



for any < 4> G C^°(Qt) and any k G R. From this we conclude easily that (j3,8[) holds. 

It remains to prove the parabolic entropy inequality ()3.9|) , Let < cp G Cq°(Qt) and 
k G R be such that A{k) ^ H. Starting off by choosing \sgn 11 (A(w £ ) — A(k))(f>\ as a test 
function in the weak formulation and then continuing exactly as in the proof of ()3.3j) , we 
obtain 

E par (w £ , k,<f>) = lim J J eVw £ ■ V[sgn v (A(w £ ) - A(k))<f>] dtdx. 
Qt 

The right-hand side of this equality can be expanded into 
lim J J (eA'(w £ )\Vw £ \ 2 sgn' v (A(w £ ) - A(k))4> 



Q 



T 



+ e sgn v (A(w £ ) - A(k))Vw £ ■ V</>J dt dx 

> lim J J esgn v (A(w £ ) - A{k))Vw £ -Vcpdtdx > -e j j \Vw £ ■ V0| dtdx. 
Qt Qt 
This concludes the proof of (|3.9|) . ■ 

4 Proof of Theorem [T] 

Following Carrillo [5] (see also ^3]) m this section we use Lemma |2] and Lemma |3] to 
prove Theorem ^ Let w £ = w £ (x,t) solve Ql.l|) and w = w(y,s) solve ()1.2(l . Following 
Kruzkov and Kuznetsov 16 a we now specify a nonnegative test function <j> = <j)(t, x, s, y) 
defined on Qt x Qt- To this end let p G Cq°(R) be a function satisfying 

supp(p) C {cr G R : |ct| < 1}, /9(<t) > Vct G R, / p(a) da = 1. 

J®. 

For x G R d , i £ R and r,r > 0, let w r (x) = \p {^) ■ ■ ■ \p {^) and p ro (t) = ±p (± 
Pick any two points u, r G (0, T), ^ < r. For any ao > define 

tp ao (i) = (t-v)- (t - t) , H ao (t) = / p ao (£) d£. 

J —oo 
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With < ro < mm(u, T — r) and ao G (0, min(f — ro, T — r — ro)) we set 

<f){x,t,y,s) := ij) ao (t)uj r (x - y)p ro (t - s). (4.1) 

Note that supp(</>(x, •, y, s)) C (ro, T — ro) for all x, y G M d , s G (0, T) and supp(^>(x, t, y, •)) 
C (0, T) for all x,y £ R d , t G (0, T). Consequently (x, t) i-> i, y, s) belongs to C^(Q T ) 
for each fixed (y, s) G Qt and (y, s) cf)(x,t,y,s) belongs to Cq 3 (Qt) for each fixed 
(x,t) G Q T . 

Observe that with the choice of 4> as in (|4.1j) we have 



d t 4> + <9 S = [p ao {t -v)- p ao (t - r)]u r (x - y)p ro (t - s), 



V x 6 + V. 



(4.2) 



Before continuing we need to introduce the two "hyperbolic" sets 

H £ = [(x,t) G Q T : A(w £ (x,t)) G h}, H = {(y,s) G Q T : A(™(y,s)) G h} 



and note that 



V x A(w £ ) = a.e. in H £ and V y A(w) = a.e. in H, 
sgn(w £ — w) = sgn(^4(ii; e ) — A(w)^j 

(Qt \ H) x Qt] U [q t x (Qt \ W £ 



a.e. m 



(4.3) 



(4.4) 



Using the approximate hyperbolic entropy inequality (|3.8j) for the viscous solution 
w £ = w £ (x, t) with k = w(y, s), we get for (y, s) G Qt 



w — w 



d t <j> + sgn(u; £ - w) [V(x)(f(w £ ) - f(w)) - V x A(w £ )] ■ V x <f> 



— sgn(w £ — w)dw x V(x)f(w)4>) dtdxdsdy > —R 



(4.5) 



Using the approximate parabolic entropy inequality (|3.9j) for the viscous solution w £ 
w £ (x, t) with k = w(y, s), we get for (y, s) G Qt \ T~t 



w — w 



d t <f> + sgn(u; £ - w) [V(x)(f(w £ ) - f(w)) - V x A(w £ )] ■ V x <f> 



Qi 



sgn(w £ — w)div x V(x)f(w)(f>) dtdx 



> lim // \V x A(w £ )\ 2 sgn'JA(w £ ) - A(w))<pdtdx - R v 



(4.6) 



Qt 



Next we would like to integrate (|4.5j) and (|4.6|) over (y, s) G Qt and (y, s) G Qt \ H. 
respectively. To this end we need to know that the involved functions are (y, s) inte- 
grable. Consider first (y, s) i— ► JJ sgn(v — u)V x A(w £ ) ■ V x (pdt dx. We denote this function 



by D(y,s). 
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To see that D(-, •) is integrable on Qt we observe that for each fixed (y, s) 6 Qt 
sgn(v - u)V x A(w £ ) = V x \ A{w e ) - A(w)\ for a.e. (x,t) € Qt 
and hence 



D(y,s) = JJ\v x \A{w £ )-A{w)\ 



• V x <fr dt dx. 



Since the function (x,t) i-> <j>(x,t,y,s) belongs to Cq°(Qt) for each fixed (y,s) 6 Qt, an 
integration by parts in x gives 



£>(y,s) = - J J \A(w e ) - A(w)\A x cf)dtdx. 

Qt 



Integration over (y, s) G Qt and estimation yield 



J J D(y,s)dsdy < JJJJ (\A{w £ (x,t))\ + \A(w(y, s))|) A x <j>(x,y,t, s) dt dx ds dy. 



By changing the variables (z := x — y, r = i — s) and taking into account that w £ ,w £ 
L x (Qt) we find that 

J J D(y,s)dsdy < JJJJ \A(w £ (x,t))\^j ao {t)\A z uj r (z)\ Pro (r) dt dx dr dz 



-+ nil \A(w(x - z,t-T))\ip ao (t)\A z uj r (z)\p rQ (T)dtdxdTdz 
< \\A(w £ )\\ L i {QT) \\A z uj r \\ L i {Rd) + \\A(w)\\ L i {QT) \\A z uj r \\ L i {Rd) < oo. 

Hence we have that D(-,-) is integrable on Qt- 

In a similar vein one can also show the integrability of 

(y, s) i— > yy |tt; £ — u> \dt4>dt dx, 
Qt 

(y,s)^ J J sgn(w £ -w)V(x)(f(w £ ) -/(«;)) • V x (pdtdx, 
Qt 

(y, s)i-> sgn(w £ - w)dw x V(x)f(w)4>dtdx, and (y, s) i-> i? visc . 



It remains to consider the integrability of the function 

Q T \ri3 (y,s)^lim J J \V x A{w £ )\ 2 sg^{A{w £ ) - A{w))<j) dt dx. 



Qi 
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This follows from (|4.fi|) . We have by Lebesgue's dominated convergence theorem and the 
first part of (|Q)l 



J I Aim J j \V x A(w £ )\ 2 sg^ v {A{w s ) - A(w))4>dtdx \ dsdy 



}t\H Qt 



lim //// \V x A(w £ )\ sgn' v (A(w £ ) - A(w))(j) dtdxdsdy. 



(Q T \H)xc 



T 



= lim JJJJ \V x A(w £ )\sgn' v {A(w £ )-A(w))(l)dtdxdsdy. (4.7) 
(Qt\H)x(Q t \H) 

We now integrate (|4.5|) over (y,s) G Qt and (|4.6() over (y,s) S Qt \ H. Addition of 
the two resulting inequalities yields 

w £ - w\d t <f> + sgn(w £ - w) [V(x)(f(w £ ) - f(w)) - V x A(w £ )] ■ V x <f> 




QtxQt 



— sgn(w £ — w)div x V(x)f(w)(f)j dtdxdsdy 

\w £ - w\d t (f) + sgn(w £ - w) [V(x)(f(w £ ) - f(w)) - V x A(w £ )] ■ V x cp 




(Q T \n)xQ T 

— sgn(w £ — w)div x V(x)f(w)(j) S j dt dx ds dy 

+ JJJJ (K - w \ d t<t> + sgn(u; £ - w) [V(x)(f(w £ ) - f(w)) - V x A(w £ )] ■ V x <k 
HxQt 

— sgn(w £ — w)drv x V(x)f(w)(p\ dt dx ds dy 

>lim JJJJ \V x A(w E )\ 2 sga' v (A(w E ) - A(w))<f> dtdxdsdy -R visc , (4.8) 

(Q T \n)x(Q T \Hc) 

where -R v isc '■= JJ Rvisc ds dy and we have used (|4.7|) . 

Qt 

Similarly, using the hyperbolic, parabolic entropy inequalities (|3.1j) . 1)3. 3 j) for the exact 
entropy solution w = w(y,s) with k = w £ (x,t) and then integrating over (x,t) £ Qt, we 
get 

JJJJ{\ W ~ W< W + S ^( w ~ \y{v)UW - fW) - VyA(w)} ■ V y <j> 
QtxQt 



sgn(w — w £ ) div y V (y) f (w £ )4>) dt dx ds dy 



>lim JJJJ \VyA{w)\ 2 sgn v {A{w) - A{w £ ))(f) dtdxdsdy. (4.9) 

(Qt\W)x{Q t \H) 
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Using (jOj) and we find that 




QtxQt 

= - JJJJ sgn(A(w e ) - A(w))V x A(w e ) ■ V y (f) dt dx ds dy 
Qtx(Qt\W) 

= - lim J I J j sgu. v (A(w e ) - A(w))V x A(w £ ) ■ V y 4> dt dx ds dy 

Q T x(Q T \W) 



lim JJJJ V y A(w) ■V x A(w £ )sgn r) (A(w £ ) - A(w))<f) dt dx ds dy. 
Qtx(Qt\H*) 




= -lim JJJJ V y A(w) ■ V x A(w £ )sgn' v (A(w £ )- A{w))<pdtdxdsdy. (4.10) 
{Qt\H)x{Q t \W) 

Similarly, again using (|4.3|) and 1)4.4(1 . we find that 

QtxQt 

= -lim JJJJ V x A(w £ ) ■ V y A(w)sgn' v (A(w)- A(w £ ))<j) dt dx ds dy. (4.11) 

(Qt\H)x(Q t \w) 

The use of the second part of (|4.2|) when adding (|4.8|) and (|4.10f) yields 

\w £ - w\d t 4> + sgn(u> £ - w) [V(x)(f(w £ ) - f(w))] ■ V x (f> 
QtxQt 

— sgn(w £ — w)div x V(x)f(w)<f) S J dt dx ds dy 

> lim JJJJ (\V x A(w £ )\ 2 - V y A(w) ■ V x A(w £ )^j 

(Qt\W)x(Q t \H) 

x sgn v (A(w £ ) - A(w))(j) dt dx ds dy - R visc . (4.12) 
Similarly the addition of (|4.9|) and (|4.11|) yields 

\w - w £ \d s (f> + sgn(w - w £ ) [V(y)(f(w) - f(w 6 )} ■ V y (f> 

QtxQt 





— sgn(w — w £ ) div y V (y) f (w £ )4>\ dtdx ds dy 



>lim JJJJ [\V y A{w)\ 2 -V x A(w £ )-V y A(w] 

(Qt\W)x(Q t \H) 

x sgn'^Aiw) - A(w £ ))(j)dtdxdsdy. (4.13) 
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Following Karlsen and Risebro we write 

sgn(u> £ — w)V(x)(f(w £ ) — f(w)) ■ V x (f> — sgn(u; e — w)div x V(x)f(w)(j) 

= S gn(w £ - w){V(x)f(w £ ) - V(y)f(w)) ■ V x </> 

+ sgn(w £ - w)div x [(V(y)f{w) - V (x) f (w)) 4>] , 
sgn(u> - w £ )V(y)(f(w) - f(w £ )) ■ V y (p - sgn(w - vf)div v V(y)f(uf)<j> 

= sgn(w £ - w)(V(x)f(w £ ) - V{y)f{w)) • V y 

- sgn{w £ - w)dw y [(V(x)f(w £ ) - V(y)f(w £ ))4>] . 

When adding (|4.12|) and (|4.13j) . we use the second part of (|4.2|) and the identities 
sgn(— r) = — sgn(r) a.e. in R, sgn^(— r) = sgn^(r) a.e. in R. 
The final result takes the form 



\w £ — w\(dt4> + d s (j>) dt dx ds dy 
QtxQt 

5; -^diss ~t~ -^visc Rconv — -^visc -^conv) (^-l^) 

where the expression for dt<ft + d s 4> is written out in (|4.2jl . 

-Rconv := JJJJ I conv dtdxdsdy, 
QtxQt 



/ conv := sgn(w £ - w)[div x [(V(y)f(w) - V (x) f (w)) </>] 
-div y [{V(x)f(w £ )-V(y)f(w £ 



and 



R diss := -Hm / / / / IV^K) - V y A(w)\ 2 




(Qt\Hc)x(Q t \H) 

x sgn„ (A(w e ) - A(w))4> dtdxdsdy < 0. 
Having in mind the first part of (|4.2j) . we get by the triangle inequality 

\w £ (x,t) - w(y,s)\(d t (p + d s (f)) dtdxdsdy < R w e jW + R W:X + R W)t , 
QtxQt 

where 

Rw*,W ■= - JJ jj \w £ (x,t) - W(x,t)\ [p ao (t -V)- p ao (t - T)] 
QtxQt 

x uj r (x — y)p ro (t — s) dtdxdsdy, 

R w ,x ■= ~ JJJJ \w(x,t) - w(y,t)\ [p ao (t -v)- p ao (t - t)] 
QtxQt 

x Lu r (x — y)p ro (t — s) dtdxdsdy, 
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Rw,t ■ = ~JJJJ \ W (y^) ~ U, (y' S )l[Pao( t - v ) - Paoit-T)] 
QtxQt 

x u r (x — y)pr (t — s) dtdxdsdy. 
Firstly a standard L 1 continuity argument gives lim R w t = 0. Next 

r |0 



lim R wx 

aolO 



w(x, t) — w{y, t)\ — \w(x, v) — w(y, v)\\oj r (x — y) dx dy 



(z:=x~y) 

< I I \w(y + z,t) - w(y,T)\oj r (z)dydz 



< \ w \L°°(0,T:BVQBL d )) / \ z \^r{z) dz < C X r, 

where C\ := | , w|£ oo (0,T;.BV(]R< i ))- Finally we have 



lim Rw e ,w = / \w £ (x, t) — w(x, r)| dx — / \w e (x, v) — w(x, v)\ dx. 
Q oJ.O J R d J R d 

In summary from <|4. 14|) we obtain the following approximation inequality 

\w £ (x, r) — w(x, t)\ dx 

< / \w £ (x, v) - w(x, v)\ dx + C\r + lim ( R visc + -R con v ) • (4-15) 

JU d r ,ao|0\ / 

We start with the estimation of R v i sc , which can be done as follows: 
d 



~ V ) | Pro — s) dtdx ds dy 

i=1 Q T xQ T 

^ AT t* f 

y2 / \dxiW e \\d Xi uJr{x -y)\dxdydt 

~1 Jv JR d JR d 



aolO 

i=l 



<eK/ r y2 / \d Xl w £ \ dtdx <eTK/r\w £ \ Laa( n T . Bvm) < C 2 Te/r, (4.16) 

4 = 1 J Rd 

where K := J Rd \5'(a) \ da and C 2 := K\w £ \ L oo^ T . BV( ^d)y 

Before we continue with the estimation of R conv we write I CO nv = -^conv + -^conv where 

/i onv = sgn(w £ - w) [{V(y)f(w) - V(x)f(w)) ■ V x </> 

- [V(x)f(w £ )-V(y)f(w £ )) - V y <j>], 
/ c 2 onv = sgn(u> £ - w)(div y V{y)f(w £ ) - div x V(x)f(w))^, 
so that 

-^conv = / / / / ^conv dt dx ds d y, R conv = J conv dt dx ds d V- 



'conv "ai -'"conv I I I I conv 

QtxQt Qt~xQi 
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We start by estimating Rl onv - To this end introduce 
F(w £ ,w) := sgn(w £ — w) [f(w 6 ) — f(w)] 



and observe that since V y (j) = —V x (p, 



«conv = //// {( y ( x ) ~ V(y))F(w £ ,w)^ ■ V x 4>dtdxdsdy. 



J t xQt 



The function F(-,-) is locally Lipschitz continuous in both variables and the common 
Lipschitz constant equals Lip(/). Since w £ G L°°(Q T ) n L°°(0, T; BV(R d )), V x F(w £ ,w) 
is a finite measure and 

d Xi F(w £ , w)\ dt dx < Lip(/) \d Xi w £ \ dt dx , i = l,...,d. 



Integration by parts thus gives 

-^conv = -J I J J (div x V(x)F(w £ ,w)ip ao (t)uj r (x - y)p ro (t - s) dt dx ds dy 

QtxQt 



(V(x) - V(y)) ■ V x F(w £ ,w)*p ao (t)uj r (x - y)p ro (t - s)dtdxdsdy . 



Jt^Qt 



1 2 

For .Rconv we calculate as follows: 



|-RcoL| < L ip(/)X] JJJJ\ V ^ X) ~ Vi ( y }\ \ d XiW e \il} ao (t)ur(x - y)pr (t - s)dtdxdsdy 
i=1 Q T xQT 

' !,L ° Lip(/)V] / / / \Vi(x) - Vi(y)\ \d x .w £ \ u; r (x - y) dx dy dt 

~l Ju JR d JR d 

' ,: l.ipi./)V / / / \V i (y + z)-V i (y)\\d yt w £ (y + z,t)\uj r (z)dzdydt 

~[ Ju JR d Jl!L d 

< Lip(F)Lip(/) V / / / \z\\d yi w £ (y + z,t)\uj r (z)dzdydt 

i=1 Ju JR d Jw 1 

< TLip(V)Up(f)\w £ \ L ^ {QTBV{Rd)) / \z\u r (z)dz 

jR d 

< rTUp(V)Up(f)\w £ \ L c C(0tT . BV{Rd)) < C 3 Tr, 

where Lip(V) := max Lip(Vj) and C3 := Lip(V)Lip(f)\w £ \ L oa( QT . B ymd\\. Note that we 

i=l,...,d ' ' \ ') 

have used the Lipschitz regularity of the velocity field V (see (|2.1j0 to get the desired 
result. 
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Regarding the term -Rc 0nv , we firstly rewrite it as 



^conv = JJJJ div x V(x)F(w £ ,w)ip ao (t)u) r (x - y)p ro (t - s)dtdxdsdy 
QtxQt 



p2,l 
^com 



+ JJJJ Sgn ( w£ ~ w }(^ v y V (y) - div^y (x)) f (w £ )ip ao (t)u r (x - y)p ro {t - s) dtdxdsdy . 



2,2 
conv 



We set 

D(x) = divV(x), 

and keep in mind that D G BV{R d ) by ()2.1|) , Now we estimate -Rconv as follows: 

|-RciL| ^ \\f( w£ )\\L^(Q T ) //// \D{y) - D(x)\ip ao (i)u r (x - y)p ro (t - s) dtdxdsdy 



QtxQt 



II/(^)IIl~(q t ) / / \D(y)-D{x)\uj r {x-y)dxdydt 

Jv JR d JR d 

{Z= = V) \\f(w £ )\\ L oo {QT) 1^ f f \D{y)-D(y + z)\u T (z)dzdydt 



< T \\f( wS )\\L°°(Q T )\D\ BV{Rd) J ^ \z\uj r (z)dz < C 4 Tr, 

where C4 := ||/(tD £ )||^oo(Q T )|L'| B ^^ Rd ^. Note that we have used the BV regularity of divV 
to get the desired result. Since Rconv — Rcon-vi we have 

Rcon, = flconv + ^conv < C 5 Tr, C 5 = max(C 3 , C 4 ). (4.17) 

Set C 6 = max(Ci,C , 2 ,C 5 ). Then from (jUSJ), (HTTCt and l|07|l we get 

/ \w £ (x, t) — w(x, r)| dx < [ \w e (x, v) — w(x, v)\ dx + C$ ( (1 + T)r H -J 

Jr* jR d \ r J 

^C 6 ((l + T)r + ^). (4.18) 
By choosing r = y/Te we immediately obtain 

/ \w £ (x,t) - w(x,t)\ dx < C 7 VTe (4.19) 

jR d 

for some constant C7 independent of e. To obtain (|2.5[) . we simply integrate (|4.19|) over 
r€(0,T). 
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Added in process 

After the main result of this paper was obtained, we became aware of a paper by Eymard, 
Gallouet and Herbin ^U] which also proves an error estimate for viscous approximate 
solutions. They, however, deal with a certain boundary value problem with a divergence 
free velocity field and obtain an error estimate of order es. As is the case herein, the proof 
in jlOj does not rely on a continuous dependence estimate. 
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